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1. Introduction 

The Darboux-Crum transformation is a well-known and powerful technique in Quantum 
Mechanics for generating new exactly solvable potentials from known ones [121 US] • A 
situation of particular interest arises when the Darboux transformation can also be 
applied to construct new families of orthogonal polynomials from known ones, since it is 
often the case that the bound states of an exactly solvable potential are polynomial 
after a change of independent variable and a rescaling of the wave function by a 
suitable non-vanishing weight function. Of course, care has to be exercised in order 
to characterize those cases in which the eigenfunctions obtained by this procedure 
do indeed give rise to orthogonal polynomial families which will in turn correspond 
to a well-defined Sturm-Liouville problem. In particular, once it is known that the 
transformed eigenfunctions are polynomial, one still has to show that they are complete 
in the underlying weighted L 2 space. Our purpose in this paper is to apply these ideas 
to generate novel families of complete orthogonal polynomial systems that are solutions 
of Sturm-Liouville problems, and in particular to show how the families of generalized 
Laguerre orthogonal polynomials recently constructed by Odake, Sasaki [23l l24l l25l 26j. 
and Quesne [271 EE] fit into the classification principle for Darboux transformations first 
introduced in [12]. It is worth to stress that these new exceptional polynomial families, 
although solutions of a Sturm-Liouville problem, are outside the Askey- Wilson class [2]. 

To set our results in context, we first recall the foundational theorem of Bochner 
[5] which states that if an infinite sequence of polynomials {P n (z)}^=o satisfies a second 
order eigenvalue equation of the form 

p(x)P' r [(x) + q{x)P' n {x) + r{x)P n {x) = X n P n (x), n = 0, 1, 2, . . . (1) 

then p(x),q(x) and r(x) must be polynomials of degree 2,1 and respectively. In 
addition, if the {P n (x)}^ =0 sequence is an orthogonal polynomial system, then it has 
to be (up to an affine transformation of z) one of the classical orthogonal polynomial 
systems of Jacobi, Laguerre or Hermite pfl 1221 El EQj [19] . 

In a pair of recent papers [TH [T7J, we have shown that there exist complete 
orthogonal polynomial systems, defined by Sturm-Liouville problems, that extend 
beyond the classical families of orthogonal polynomials arising from Bochner's 
classical theorem on the characterization of Sturm-Liouville polynomial systems. 
What distinguishes our hypotheses from those made by Bochner is that the first 
eigenpolynomial of the sequence need not be of degree zero, even though the full set of 
eigenfunctions still forms a basis of the weighted L 2 space. The situation we considered 
in [TTJ is that of complete orthogonal polynomial systems starting in degree one. For 
this m = 1 case a full characterization of all Sturm-Liouville polynomial systems 
is available thanks to the classification of codimension one exceptional polynomial 
subspaces performed in [16]. The concept of an exceptional polynomial subspace 
was introduced in [HI [T5] . For some recent applications of exceptional orthogonal 
polynomials see [211 EI] • 
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In the present paper, we pursue this program in higher codimension by constructing 
explicit examples of complete systems of higher codimension using the Darboux 
transformation in a systematic fashion. In particular, we will construct the analogues in 
codimension m of the Xi Laguerre polynomials. Our paper is organized as follows. In 
Section 2, we define the notion of an m-orthogonal polynomial system arising from 
modules of higher codimension and we introduce the notion of polynomial Sturm- 
Liouville problem, corresponding to the case in which the eigenfunctions of the Sturm- 
Liouville operator are polynomials. In Section 3, we carefully study the various 
cases in which the Darboux transformation preserves the polynomial character of the 
eigenfunctions (these are precisely the algebraic Darboux transformations introduced 
in [E]). We also show the role played by the property of shape invariance in the 
various factorizations that give rise to Darboux-Crum transformations. In Section 4, 
we apply these results to the case of the Sturm-Liouville problem defining the Laguerre 
polynomials and show precisely how the LI and L2 families of codimension m Laguerre 
polynomials obtained in (25] fit into our general classification scheme. In Section 5, we 
show that these polynomials satisfy remarkable shape invariance properties that arise 
from the intertwining relations obtained between the partner second-order operators 
and iterations of the first-order operators defining the Darboux transformations. We 
will see that the shape invariance reported in [23] [27] follows from the shape invariance 
of the initial operator and the formal properties of the Darboux-Crum transformations. 

The perspective taken in this paper is largely that of the formal calculus of 
differential operators. Previously, analytic aspects of the Darboux-Crum method for 
Sturm-Liouville systems were considered in [10]. Our approach is different in that we 
focus on algebraic properties and exact solutions. In a subsequent paper [18] we shall 
study the codimension m Laguerre polynomial system in a functional analytic setting by 
giving a spectral theoretic characterization of the codimension m Laguerre polynomial 
system in the context of Sturm-Liouville theory. A detailed analysis of the asymptotic 
properties of the zeros of these polynomials will be given. We shall also further develop 
some of the key formal properties of these polynomials that result from the factorization 
and shape invariance of their defining operators, namely their orthogonality properties, 
Rodrigues-type formulas and generating functions. Finally we mention that this entire 
analysis can also be carried out in the case of Jacobi polynomials. 

2. Preliminaries 

We will say that a differential operator 



is exactly solvable by polynomials (PES) if it admits infinitely many real, polynomial 
eigenfunctions Uj(x): 




(2) 



T (Vj) = ^jVj, degyj < degy j+1 , Xj G R, j = 1, 2, . . . . 



(3) 
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Moreover, we say that a sequence of polynomials {yj}fLi has codimension m if 

deg yj = m + j - 1, j = 1, 2, . . . (4) 

We also say that T(y) is an m-PES operator if the eigenpolynomials satisfy the above 
condition. We say that T(y) is primitive, if the eigenpolynomials do not possess a 
common root (real or complex). 

Note: if T(y) has at least 3 linearly independent polynomial eigenfunctions, then 
necessarily, p(x),q(x),r(x) must be rational functions. 

Let / = (xi,X2) be an open interval (bounded, unbounded, or semi-bounded) and 
let Wdx be a positive measure on / with finite moments of all orders. We say that a 
sequence of real polynomials {%}^Li forms an orthogonal polynomial system (OPS for 
short) if the polynomials constitute an orthogonal basis of the Hilbert space L 2 (J, Wdx). 
If @ holds, we speak of an m-OPS. 

The following definition encapsulates the notion of a system of orthogonal 
polynomials defined by a second-order differential equation. Consider a boundary value 
problem 

- (PyJ + Ry = XWy (5) 
lim (Py'u - Pu'y)(x) = 0, z = 1, 2, (6) 

x-+xf 

where P(x),W(x) > on the interval / = (xi,x 2 ), and where u(x) is a fixed 
polynomial solution of (jSJ). We speak of a polynomial Sturm-Liouville problem (PSLP) 
if the resulting spectral problem is self-adjoint, pure-point and if all eigenfunctions are 
polynomial. We speak of an m-PSLP if the eigenpolynomials satisfy (T4]). If m = 0, then 
we recover the classical orthogonal polynomials, the totality of which is delineated by 
Bochner's theorem. For m > 0, Bochner's theorem no longer applies and we encounter 
a generalized class of polynomials; we name these exceptional, or X m polynomials. 
Given a PSLP, the operator 

T{y) = W-\Py')' -W- x Ry 

is PES. Letting p(x),q(x),r(x) be the rational coefficients of T(y) as in (j2|), we have 



P{x) = exp q/pj , (7) 

W(x) = {P/p)(x), (8) 
R(x) = -(rW)(x), (9) 

Hence, for a PSLP, P(x), R(x), W(x) belong to the quasi-rational class [TTj. meaning 
that their logarithmic derivative is a rational function. 

Conversely, given a PES operator T(y) and an interval I = (x\,X2) we formulate 
a PLSP (jSJ) by employing dTJ)-© as definitions, and by adjoining the following 
assumptions: 

(i) P(x), W(x) are continuous and positive on / 

(ii) Wdx has finite moments, i.e. L x n W(x)dx < oo, n — 0, 1,2, . . . 
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(iii) \im x ^ Xi P{x)x n = 0, i = l,2, n = 0,l,2... 

(iv) the eigenpolynomials of T(?/) are dense in the Hilbert space L 2 (7, Wdx). 
These definitions and assumptions (i) and (ii) imply Green's formula: 

/ T(f)g Wdx - / T(g)f Wdx = P(fg - fg') (10) 

By (iii) if f(x),g(x) are polynomials, then the right-hand side is zero. If / and g are 
eigenpolynomials of T(y) with unequal eigenvalues, then necessarily, they are orthogonal 
in L 2 (J, Wdx). Finally, by (iv) the eigenpolynomials of T(y) are complete in L 2 (7, Wdx), 
and hence satisfy the definition of an OPS. 

We now describe a construction that systematically generates polynomial Sturm- 
Liouville systems of arbitrarily large codimension m > 0. 

3. The Darboux transformation 

Let T(y) be a differential operator (T5]) with rational coefficients. We speak of a rational 
factorization if 

T-\ Q = BA (11) 

where A(y),B(y) are first order operators with rational coefficients and where Ao is a 
constant. Let us write 

A(y)=b(y'-wy), (12) 
B(y)=b(y'-wy), (13) 

where w(x),w(x), b(x), b(x) are all rational functions. Given a rational factorization we 
introduce the partner operator 

f = AB + A (14) 

whose explicit form is 

T(y) = py" + qy' + fy (15) 

where 

q = q + p ' -2pb'/b (16) 
f = —p(w' + w 2 ) — qw + Ao- (17) 
We will refer to 



<p{x) = exp yj w 

as a factorization eigenfunction (quasi-rational) and to b(x) as the factorization gauge 
(rational). The former satisfies 

T(0) = A o 0. (18) 

Equivalent ly, 

w(x) = <f>'(x)/</)(x) (19) 
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is a rational solution of the following Ricatti-like equation: 

p(w' + w 2 ) + qw + r = A . (20) 

For a fixed T(y), a rational factorization is fully determined by a quasi-rational 
factorization eigenfunction and a rational factorization gauge. Indeed, given p, q, r, w, b 
relation (fTTj) gives us 

b = p/b, (21) 
w = -w - q/p + b'/b. (22) 

The choice of b(x) determines the gauge of the partner operator. Consider two 
factorization gauges 6i(x),6 2 (a ; ) and let T 1 (y),T 2 (y) be the corresponding partner 
operators. Then, 

T 2 = /U -1 Ti/i, 

where 

/i(x) = b 1 (x)/b 2 {x). 

The above construction of the partner operator is symmetric with respect to the 
interchange of the hatted and unhatted variables. Letting P(x) be as in (GO), and setting 

^ (x) = exp (7 = (p/b)l)<p(xy (23) 

we have 

f (0) = Xol 

We also have 

P/b = P/b, (24) 

bb = p, (25) 

q/p - b'/b = q/p - b'/b (26) 

Thus, starting with T and taking as the factorization function and b(x) as the 
factorization gauge, we recover T(y). 

Next, suppose that T is a PES operator with eigenpolynomials {%}. If /x(x) is a 
polynomial, then \xT\i~ x is also a PES operator, with eigenpolynomials {f^Vj}- Therefore, 
we can fix the gauge of a PES operator by requiring that the eigenpolynomials are 
primitive (no common roots). 

By construction, partner operators obey the following intertwining relations: 

TA = AT, BT = TB. (27) 

Hence, if T is a PES operator with eigenpolynomials {%•}, then {A(yj)} are 
eigenf unctions of the partner operator T with the same eigenvalues. By inspection 
of ([121) . with the appropriate choice of b(x), the A(yj) are polynomials. Hence, if T 
is PES, then so is T. Furthermore, the requirement that the eigenpolynomials of T 
be primitive fixes b(x) up to a choice of scalar multiple. In many cases, such as the 
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factorization shown in (157|) - (l69p . it will suffice to take b(x) to be the denominator of 
w(x). However, there are other cases, such as the factorization shown in (I128l) - (ll3ip . 
where b(x) must be a rational function. 

The duality between T and T has another aspect. Let W(x) be as in (jSj) and let 
W(x) be analogously defined. Hence, by equations (121]) ( 1251) . 

W = P/b 2 = pW/b 2 (28) 

Consequently, A and —B are formally adjoint relative to these measures: 

A(f)g Wdx + / B(g)f Wdx = (P/b)fg (29) 

XI J X\ Xl 

If the above RHS vanishes for polynomial f,g, then A and —B, with suitably defined 
domains, give rise to adjoint operators in the rigorous sense of densely defined linear 
operators on Hilbert spaces L 2 (J, Wdx) and L 2 (J, Wdx), respectively. 

Darboux transformations can be classified into three types as far as their spectral 
properties are concerned [El [30]: state-deleting, state-adding, or isospectral. 

(i) state-deleting transformation: In this case the factorizing function <f>(x) satisfies 
G L 2 (I,Wdx) and the formal factorizing eigenvalue Ao is the maximu of the 
spectrum of T. 

(ii) state-adding transformation: In this case the partner factorizing function <p(x) 
satisfies <j) G L 2 (I,Wdx) and the formal factorizing eigenvalue Ao must be above 
the maximum of the spectrum of T. Equivalently, from ( 1231) and ( 1281) it follows 
that 

<p G L 2 (I, Wdx) -y^ 1 e L 2 (I, Wdx) 

so it is clear that the spectral properties of the transformation only depend on the 
choice of 0, not on the choice of gauge b(x). 

(hi) isospectral transformation: In this case (j) ^ L 2 (I ,Wdx), <fi ^ L 2 (I, Wdx) and 
the formal factorizing eigenvalue Ao must be above the maximum of the spectrum 
of T. 

In the context of algebraic Darboux transformations discussed in this paper, if 
we assume that both T and T are PSLPs, the above spectral characterization can be 
particularized to a purely algebraic one. 

(i) A state-deleting transformation corresponds to <p — yi, the first eigenpolynomial of 
T. 

(ii) A state- adding transformation corresponds to (ft (as defined by (1231 ) being a 
polynomial. 

| Note that, as opposed to the usual convention in Schrodinger operators where the spectrum is 
bounded from below, in this paper the spectrum of all Sturm-Liouville problems is bounded from 
above. The eigenfunction corresponding to the maximum of the spectrum corresponds therefore to the 
ground state. 
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(iii) Isospectral transformations correspond to neither nor <f> being polynomials. 

The above conditions can be explicitly verified on the factorizations performed in 
Sections 4 and 5. For example, equations (!67|) -( l69|) show an isospectral factorization; 
neither of the factorizing eigenf unctions is a polynomial. By contrast, equations 
(1 1 2 81) - ( 1X311) show a state-deleting/state-adding factorization; one of the factorizing 
eigenfunctions is a polynomial, while its partner eigenfunction is not. 

State-adding and state-deleting factorizations are dual notions, in the sense that if 
the factorization of T is state-deleting, then the factorization of T is state-adding, and 
vice versa. 

As we already pointed out, the eigenpolynomials {yj} and {yj} constitute 
orthogonal polynomial systems relative to L 2 (7, Wdx) and L 2 (7, Wdx), respectively. 
The adjoint relation between A and B allows us to compare the L 2 norms of the two 
families. Indeed, by (TTT|) (13|) (|29|), 

f{A{ yj )f Wdx = -jB(A( yj )) yj Wdx = (A - A,) Jy 2 3 Wdx (30) 

3.1. Shape-invariance 
Suppose that 

T k (y) = p(x)y" + q k (x)y' + r k (x)y, k e K, (31) 

is a family of PES operators, where K is some parameter index set. If this family is 
closed with respect to the state-deleting Darboux transformation, we speak of shape- 
invariant operators and polynomials. To be more precise, let Ti k (x) = yk,i(x) be be 
the corresponding ground-state eigenpolynomial. Without loss of generality, we assume 
that that the ground-state energy is zero, and let 

T k = B k A k , A k {ix k ) = (32) 

be the corresponding factorization. Shape-invariance means that there exists a one-to- 
one map h : K — >■ K and real constants such that 

Th( k ) = A k B k + X k . (33) 
Necessarily, there exist constants a k j,f3 kt j such that 

Vh(k),j-i = ac ktj A k (y ktj ), j > 1, (34) 

Vkj+i = 0k,j+iB k (y h (k),j), j > 0, (35) 

Pk,jOc k:j = X kJ (36) 

^h( k )j = Afcj+i + X k . (37) 
In accordance with (TTI), define 



P k (x) =exp (J q k /pj (38) 

Let b k (x) denote the shape-invariant factorization gauge; i.e.; 

A k {y) = {b k /7r k )W{TT k ,y), (39) 
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W(/, g) = fg' - f'g (40) 

Equation ( {TBI implies the following necessary condition, 

P Pk/P m = b 2 k . (41) 

This is a rather strong constraint, because the left-hand side is a product of quasi- 
rational functions, while the right-hand side is a rational squared. 

3.2. Covariant factorization 

Next, we introduce the notion of a covariant isospectral factorization. Let T k (y) be a 
shape- invariant family of PES operators, as above. Suppose that <p k (x) is an indexed 
family of isospectral factorization functions. Let 

T k = B k A k + X k , A k {<f) k ) = (42) 

be the corresponding isospectral factorization. Let 

f k = A k B k + A fc , B k (k) = (43) 

be the partner operator and partner eigenfunction, respectively. We say that the 
factorization with respect to <p k is covariant if 

A k (cj) k ) oc (j) h{k) . (44) 

The following Lemma furnishes a useful test for covariant factorization. Let us say that 
a PES operator is formally non-degenerate if for every formal eigenvalue (eigenfunction 
is quasi-rational) Ao G M there exists at most one linearly independent quasi-rational 
eigenfunction with that eigenvalue. 

Lemma 3.1 Suppose that <j) k (x) is continuous with respect to k and formally non- 
degenerate for generic values of k G K . Furthermore, suppose that 

^h(k) = Afc + Afc. (45) 

Then, the factorization with respect to <p k is covariant. 

Proof. By (J32]) (J33D and g5l), 

T m {A k {<f) k )) = A k {T k {<f) k )) + X k A k {<p k ) (46) 
= (Afc + Afc)Afc(0fc) (47) 
= \ m A k ((f) k ) (48) 

Since <f> k is quasi-rational and since A k has rational coefficients, A k (<p k ) is also quasi- 
rational. Hence, (j4"4"l) holds for generic k. Therefore, it holds for all k. QED 
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4. Laguerre polynomials 

Let us introduce the PES operator 

C k (y) := xy" + (k + 1 - x)y' . (49) 

The classical associated Laguerre polynomials, L k>n (x) can be defined as the 
corresponding eigenpolynomials, 

C k (L k , n ) = -nL kiTl , (50) 

normalized by the condition 

(-l) n 

L kn (x) = ; — x n + lower degree terms. 

n! 

The classical Laguerre polynomials are shape-invariant by virtue of the following 
factorizations 

C k = B k A k (51) 
C k+ i = A k B k + 1 where (52) 

Mv) = v' (53) 

B k (y) = xy' + (k + 1 - x)y (54) 
For k > — 1, the resulting polynomials are orthogonal relative to the weight 

W k (x) = x k e~ x , xe(0,oo), (55) 

and can be realized as solutions of a spectral problem (3j [7]. The corresponding L 2 
norms are given by 

L 2 kn W k dx = Y{n + k + l)/n\. (56) 

The quasi-rational eigenfunctions of C k (y) are known [U Sec. 6.1]: 

0i(x) = L k>m (x), A = -m (57) 

<j) 2 (x) = x~ k L_ k ^ m (x) A = k - m (58) 

3 (x) = e x L k , m (-x) X = k + l + m (59) 

04 (x) = x~ k e x L_ ktm (—x), A = m + 1, (60) 

where m — 0,1,2, . . .. The corresponding factorizations were analyzed in [12]. Of these, 
0! with m = corresponds to a state-deleting transformation and underlies the shape- 
invariance of the classical Laguerre polynomials. For m > 0, the <j>\ eigenfunctions 
yield singular operators and hence do not yield novel orthogonal polynomials. The 04 
family results in a state-adding transformation. The resulting orthogonal polynomials 
do not satisfy condition (@J; such factorizations were discussed in [T3]. The type 2 and 
type 3 factorizations 02, 03 result in novel orthogonal polynomials, although for 3 it is 
necessary to assume that k > m. These families correspond, respectively, to the type 
LI, L2 Laguerre polynomials of [25] . 

Let us consider these two families of factorization on a case-by-case basis. The 
derivations that follow depend in an elementary fashion on the following well-known 



o 



X m -Laguerre polynomials 11 

identities of the Laguerre polynomials. We will apply them below without further 
comment. 

L k , (x) = 1 (61) 
L M (x) = 0, n < -1, (62) 
nL k , n (x) + (x - 2n - k + l)L Kn _i(x) + (n + k - l)L Kn _ 2 {x) = 0, (63) 
L 'k,n( x ) = -£fc+i,n-iO), (64) 

Lk,n\ x ) — Lk+l, n (x) — ^+l,n-l(4 (65) 

^.1. T/ie LI family 

Fix an integer m > and a real fc > — 1. Take 03 (x) as the factorization function and 
€k,m{x) = L k , m {-x) (66) 

as the factorization gauge. Applying (fl2|) (Fl3|) (I2TT) ([521 . the resulting factorization is 
£fc = m A km + k + m + 1, where (67) 

A k,m(y) = £k,mV' ~ €k+l,mV ( 68 ) 

Bl m (y) = (xy' + (l + k)y)/^ m (69) 

The partner eigenfunction is <p{z) = z~ l ~ k . Let us define 

£fc,m = ^k-l,m^k-l,m + k + m (70) 
£k,m(y) = x v" + ( k + 1 ~ + m y ~ 2Pk-i,m(xy' + ky), where (71) 

Pk,m £,k,m/ '^k,m £fe+l,m— l/£fe,m- (^2) 

On the basis of the above factorization, we define type I exceptional Laguerre 
polynomials to be 

^k,m,n = ~A 1 k _ l rn (Z/fc-l,n- m ) (73) 

(74) 

By construction, these polynomials satisfy 

4, m ( L fc, m ,n) = {m- n)L\ m n , n>m. (75) 
By (flDl) and fT28|) the sequence {L| mn }^ =m constitutes an m-OPS relative to the weight 

WlJ?) = x k e- x /e k -i,m, x G (0, oo) (76) 
Using ( 13"Uj) and (15B"]) . we obtain 

(4, m ,n) 2 < m dx = (k + n)T(k + n- m)/(n - m)\ (77) 
For m = the above definitions reduce to their classical counterparts; to wit, 

L \,0,n = L k,n, (79) 

WUx) = xV (80) 
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4.2. The L2 family 

Fix an integer m > and k > m, and take 02 {x) as the factorization function. Set 

r]k,m{x) = L- k>m {x)\ (81) 

take x ?7fc jm as the factorization gauge. The resulting factorization is 

C k = B k ] m A l k ] m + (k-m), where (82) 
A k, m (y) = %Vk,my' + (k- m)r] k+hm y (83) 

B l k ] m (y) = (y' - y)hk,m (84) 

The partner eigenfunction is <p{z) = e z . Based on this factorization, we define 

£fc,m = ^fe+l.m-^fc+l.m + (k + 1 — m) (85) 

£k,m(v) =xy" + (k + l-x)y' -my + 2xa k+lim (y' -y), where (86) 

0~k,m — ~Vk,m/ r lk,m = Vk-l,m-l/Vk,m (87) 

We now define the type II X m Laguerre polynomials to be 

L k ,m,n = _ Afc+l,m (£ft+l,n-m) , Tl>m (88) 

= %Vk+l,m L k+ 2,n-m-l + (m — k — l)%+2,m -^fc+l,n-m (89) 

By construction, these polynomial satisfy 

^k,m( L l,m,n) = ( m ~ n ) L k,m,n, W > 771. (90) 

Thus, the sequence {L l k l m n }^ m constitutes an m-OPS relative to the weight 

W k n m (x) = x k e- x / V t +lim , x e (0, oo) (91) 
Using ( 156|) . we also have 

r{L\ l f W* m dx = (1 + fc + "T, 2m) r(2 + k + n - m) (92) 
Jo (n-m)\ 

As above, for m = the above definitions reduce to their classical counterparts, albeit 
the polynomials have a different normalization: 

Llo, n = -(k + l + n)L k , n . (93) 

The proof that the sets {L\ m n }™ =m and {L\ l m n }^L m span dense subspaces of the 
Hilbert spaces L 2 ([0, oo), W\ m dx) and L 2 ([0, oo), W k u m dx) will be given in a forthcoming 
publication [18J. 

5. Shape-invariance of the exceptional polynomials 

In this section we prove that above defined X m polynomials are shape-invariant. The 
explanation for this remarkable fact is the commutativity/permutability of iterated 
Darboux transformations, also known as the Darboux-Crum transformation. 

Let T (y) = T(y) be a given PES operator, and let <f>i(x), . . . , (f> n (x) be quasi- 
rational eigenfunctions. Let Ti(y) = T(y) be the partner PES operator corresponding 
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to 4>i. Now Ai(4>2) is a quasi-rational eigenfunction for Ti(y). Let T 2 (y) = T\(y) be the 
corresponding partner operator. Continue in like fashion. We arrive at the following 
chain of factorizations: 

T = B 1 A 1 + Ai (94) 

'/.; ■\ j r> J I \ r j = l,...,n-l, (95) 

= B j+1 A j+1 + A i+1 (96) 

T n = A n B n + An, (97) 

where 

(A j ---A 2 A 1 )(<j )j ) = 1 j = l,2,...,n. (98) 

In the end, we obtain the following intertwining relations: 

T B = BT n , where B = B 1 B 2 ---B n (99) 

AT = T n A, where A = A n --- A 2 A X . (100) 

By construction, 

Atf j ) = 0, j = l,2,...,n. (101) 

Hence, 

A(y) = b{x)W((f>i, ...An, v)/W(<l>i, ■■■An), (102) 

where b(x) is the higher-order rational factorization gauge, and where W denotes the 
Wronskian operator. As before, b(x) is uniquely determined (up to scalar multiple) 
by the requirement that the eigenfunctions of T n (y) constitute a primitive sequence of 
polynomials. 

The key observation is that up to sign, the above definition of T n is independent 
of the order of the factorization functions. Let us exploit this commutativity to prove 
that the above-defined X m polynomials are shape-invariant. To do so, requires that we 
consider a certain 2-step factorization. 

Let Tk(y) be a family of shape- invariant PES operators as per (I3TI) . Let ^(x) = 
Vk,i( x ) denote the corresponding ground-state eigenpolynomials, and let (1321 ( 1331) be 
the corresponding factorizations, where without loss of generality the factorization 
eigenvalue is set to zero. 

Next, let <f>k{x) be a quasi-rational eigenfunction that corresponds to a covariant, 
isospectral factorization as per (|44j) . Let T^{y) be the corresponding family of isospectral 
operators as per (l4~3l) . We claim that this family is also shape-invariant. Let 

f k = B k A k , A k (TT k ) = (103) 
be the ground-state factorization of the partner operator, where 

n k = A k {7i k ) (104) 
is the new ground-state polynomial. Our claim is that 

f m = A k B k + A fc , (105) 
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For convenience, let us set 

f k = A k B k (106) 

The 2nd order intertwining relation is 

AkT k = f k A k where (107) 

Mv) = KtfWfa, <f> k , y)/W{n k , <p k ), (108) 

and where b(x) is the 2nd order rational factorization gauge whose form is not relevant 
to our argument. There are two ways to factorize Ak, the 2nd order intertwiner: 

Ah = A k A k (109) 

A k = A h{k) A k (110) 
The 2nd equation is true because by (JUl) we have 

A h {k){M4>k)) oc A h(k) ((j) m ) = 0. (Ill) 

Hence, 

AkT k = Ah(k)Af.Tf, (112) 

= A h ( k )(Th(k) — X k )A k (113) 

= (Tfc(fc) - x k)M ( 114 ) 
Hence, by Equation (11071) . 

fkAk = (Th(k) - Afc)A- (115) 

The ring of differential operators with rational coefficients has no zero divisors. 
Therefore, the desired relation (11051) follows. 

Next, let us illustrate the above result by explicitly showing the shape-invariant 
factorization for the type I exceptional Laguerre polynomials defined in the preceding 
section. The index set consists of real k > — 1. Let us set 

n k (x) = 1, (116) 

h(k) = k + l (117) 
X k = 1 (118) 
T k (y) = C k (y) (119) 

The classical Laguerre polynomials are shape- invariant; relations (|32|) (|33|) hold, as per 
( I51j) - (l54j) . Let us fix an integer m > and set 

fc (x) = e x £ k ,m, (120) 

Mv) = A im(y) (i2i) 

B k (y) = Bl m (y) (122) 

f k (y) = 4+1,™ (123) 
~\ k = k + m (124) 

These definitions realize a particular instance of the isospectral factorizations shown 
in (1421) (]43l) . By inspection of (1571) -(1601). the operator Ck is formally non-degenerate 



X m -Laguerre polynomials 



15 



for generic k. Equation (145 p is satisfied, and hence by Lemma 13.11 the isospectral 
factorization with respect to (11201) is covariant. Therefore, the operators C\ m are shape- 
invariant. 

Next, we explicitly describe the ground-state factorization for C\ m and verify the 
shape- invariance property. To determine an explicit form for A k ,B k we make use of 
formula (j4"T|) . Here, 

q k (x) = (k + 2 - x) - 2xp k+ltm (125) 
P k (x) = exp ( [ X q k /p] = e-*x k+2 /e k , m (126) 



pp k /Pk + i = e k+1 ,je k , m (127) 

In this way, we arrive at the shape-invariant factorization 

"£fc,m = 4,m^):,mi ^-k,rn(^k,rn) = 0, (128) 

4+l,m = 4rA™ + l, ^(e^-^fc-l.m) = (129) 

where 

A*,m(v) = {ik,m/ £k-l,m){y' ~ Pk,mV) (130) 

B\, m {y) = fe-i.m/Cfe.mXV + (1 + k) y) - xy (131) 
Thus, the type I polynomials obey the following lowering and raising relations: 

Ac,m(-kfc,m,n) = ~^ J k+l,m,n-ll U > 171. (132) 

Bl,m( L k+i,m,n) = (ti + 1 - m)L l kmn+1 , n>m. (133) 
In a similar fashion, we derive the following shape-invariant factorization for the 
type II polynomials. This time, we let 

<p k (x) = x~ k r) k , m , (134) 

Afc = k — m, (135) 

f k (y) = Cl 1 ^, (136) 

A k (y) = Al m , (137) 

B k {y) = B k ] m , (138) 

q k (x) = (k - x) + 2xa Km (139) 

P k (x) = exp ( / q k /p] = e- x x k /r,l m (140) 



pPk/Pk+i = vl+i, m /vl m (141) 
Applying the formulas of section 3, we obtain the following shape-invariant factorization: 

£fc,m = P l k \rn^ k l ,m-> ^k,m (Vk+2,m) = 0, (142) 

4+l, m = KmBlm + 1, 5gm(^" 1_ W,m) = (143) 



where 



A k,m(y) = (Vk+2,m/Vk+l,m)(y' + O k +2,m V) (144) 

B k,m(y) = (Vk+i, m /Vk+2, m )(xy' + (1 + k - x) y) + (145) 

+ (Vk,m-l/Vk+2,m)xy (146) 
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The type II polynomials obey the same lowering and raising relations as in (I132p . 
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